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The magnetic ground state phase diagram of the disordered Hubbard model at half-filling is 
computed in dynamical mean-field theory supplemented with the spin resolved, typical local den- 
sity of states. The competition between many-body correlations and disorder is found to stabilize 
paramagnetic and antiferromagnetic metallic phases at weak interactions. Strong disorder leads to 
Anderson localization of the electrons and suppresses the antiferromagnetic long-range order. Slater 
and Heisenberg antiferromagnets respond characteristically difi'erent to disorder. The results can be 
tested with cold fermionic atoms loaded into optical lattices. 
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Interacting quantum many-partiele systems with dis- 
order pose fundamental challenges for theory and experi- 
ment not only in condensed matter physics [l|, d, S 0, 0] , 
but most recently also in the field of cold atoms in opti- 
cal lattices d, 0, [1, H, Indeed, ultracold gases have 
quickly developed into a fascinating new laboratory for 
quantum many-body physics; see e.g. ll|, [iSl, 
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15l . Il6l . Il7| . A major advantage of cold atoms in op- 



tical lattices is the high degree of controllability of the 
interaction and the disorder strength, thereby allowing a 
detailed verification of theoretical predictions. 

In particular, these quantum many-body systems will 
allow for the first experimental investigation of the simul- 
taneous presence of strong interactions and strong disor- 
der. This very interesting parameter regime is not eas- 
ily accessible in correlated electron materials. Namely, 
at or close to half-filling where interaction effects be- 
come particularly pronounced, strong disorder implies 
fluctuations (e.g., of local energies) of the order of the 
band-width, which usually leads to structural instabil- 
ities. These limitations are absent in the case of cold 
atoms in optical lattices where disorder can be tuned to 
become arbitrarily strong without destroying the exper- 
imental setup. Since at half filling and in the absence 
of frustration effects interacting fermions order antifer- 
romagnetically, several basic questions arise: (i) How is 
a non-interacting, Anderson localized system at half fill- 
ing affected by a local interaction between the particles? 
(ii) How does an antiferromagnetic insulator at half fill- 
ing respond to disorder which in the absence of interac- 
tions would lead to an Anderson localized state? (iii) 
Do Slater and Heisenberg antiferromagnets behave dif- 
ferently in the presence of disorder? In this Letter we 
provide answers to the above questions by calculating 
the zero temperature, magnetic phase diagram of the dis- 
ordered Hubbard model at half filling using Dynamical 
Mean-Field Theory (DMFT) [UHlifll with a geometric 



average over the disorder [21j, [22, |23|, |2J, l25| and allowing 
for a spin-dependence of the density of states (DOS). 

Antiferromagnetic (AF) long-range order is a generic 
property of interacting lattice fermions with particle-hole 
symmetry, as exemplified by the Hubbard model at half 
fillin g w ith nearest-neighbor hopping on a bipartite lat- 
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28l |. Such an instability is also highly rele- 
vant for current and future experiments in optical lattices 
[2^ , where the magnetic super-exchange energy scale has 
recently been observed in a two-component bosonic mix- 
ture [30|. The influence of disorder, e.g., due to fluctu- 
ating local potentials, on interacting quantum particles 
is subtle and leads to a remarkably rich phase diagram 
which was studied by a variety of numerical techniques 

While previous investiga- 
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tions yielded important insights into the properties of 
disordered Hubbard antiferromagnets in various regions 
of parameter space, a comprehensive study, where effects 
due to Anderson localization, genuine many-body corre- 
lations and AF order are treated within the same non- 
pcrturbative theoretical framework, did not yet exist. To 
this end we here employ the DMFT — a non-perturbative 
approach to correlated lattice fermions which accounts 
for the Mott-Hubbard metal-insulator transition (MIT) 
and magnetic ordering — in combination with a disorder 
average which is able to detect Anderson localization on 
the one-particle level [21 , 2^, 2^ . Namely, by employing 
the geometric rather than the arithmetic average over the 
disorder it is possible to determine the typical local DOS 
[3^ as a dynamical mean field within the DMFT. This 
approach was recently employed to calculate the param- 
agnetic phase diagram of the disordered Hubbard model 
24 1 and Falicov-Kimball model [2^ . Thereby it was pos- 



sible to determine the MIT due to disorder (Anderson lo- 
calization) and interactions (Mott-Hubbard transition), 
respectively, as well as the transition scenario caused by 
their simultaneous presence, within a unified framework. 
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In the absence of frustration effects the Mott-Hubbard 
MIT is completely hidden by AF long-range order [l^, 
28| . To capture this feature it is necessary to gener- 
alize the investigation and include AF solutions of the 
Hubbard model with local disorder (Anderson-Hubbard 
model), whose Hamiltonian is given by 



5^ ^ij'Aa^'Jo 



(1) 

Here tij is the amplitude for hopping between the sites 
i and j, U is the on-site repulsion, riicr ~ aL'^icr the 
local fermion number operator with (oL) as the an- 
nihilation (creation) operator of a fermion with spin cr, 
and 6; are random on-site energies. In the following wc 
work with a continuous probability distribution function 
for ei, i.e., Vie,) = e(A/2 - |ei|)/A, with 9 as the step 
function. The parameter A is a measure of the disor- 
der strength. We consider a bipartite lattice with equal 
number of fermions and lattice sites (half- filled case). In 
the absence of disorder the Hamiltonian is then explicitly 
particle-hole symmetric. 

The Anderson-Hubbard model H]) is solved within 
DMFT by mapping it onto sing le-impurity Anderson 
Hamiltonians with different e, [Tglls^l. For each random 
on-site energy e;, where i belongs to one of the sublat- 
tices s —A or B, we calculate the local Green function 
G'o-s(w,ei). From this quantity we obtain the geometri- 
cally averaged local DOS p|°°™(w) = exp [(lnpcrs(w, e^))] 
[2l|,|2i|, where yOas(w, Ci) = -ImGcrs(w, eO/^^: and (O) = 
/ deiP{ei)0{€i) denotes the arithmetic average of 0(ei). 
For comparison, the arithmetically averaged local DOS 
p™*'^(<jj) = {pas{uj,ei)) is also computed. The aver- 
aged local Green function is then obtained from the 
Hilbert transform G"g{uj) = J duj' p"g{uj') / {oj—u'), where 
a = geom (arith) denotes the geometric (arithmetic) av- 
erage. The local self-energy T,"g{uj) is determined from 
the k-integrated Dyson equation S"^(w) = w — ri"^{uj) — 
1/G"g{uj) where ri"^{uj) is the hybridization function of 
the effective Anderson Hamiltonian. The latter quan- 
tity provides the position and the resonant broadening 
of single-site quantum levels and may be interpreted as 
a molecular mean-field which describes the effect of all 
other sites within the DMFT. The self-consistent DMFT 
equations are closed by the Hilbert transform of the 
Green function on a bipartite lattice 
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where No{e) is the non- interacting DOS and s denotes 
the sublattice opposite to s. In the following we choose a 
model DOS, iVo(e) = 2VD^ - e^/nD^, with bandwidth 
W = 2D, and set VK = 1. For this DOS and a bipartite 
lattice the local Green function and the hybridization 
function are connected by the simple algebraic relation 
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FIG. 1; Magnetic ground state phase diagram of the 
Anderson-Hubbard model at lialf-filling as calculated by 
DMFT with a spin resolved local DOS (see text). 



?7^,(w) = D'^G'^s{^)/A [ia. The DMFT equations are 
solved at zero temperature by the numerical renormal- 
ization group technique [s^ , which allows us to calculate 
the geometric or arithmetic average of the local DOS in 
each iteration loop. 

To characterize the ground state of the Hamiltonian 
([T]) the following quantities are computed: The local DOS 
p'^si^) foi' ^ given sublattice s and spin direction cr, the 
total DOS for a given sublattice s at the Fermi level 
Nf{Q) = ~ '^)' staggered magneti- 

zation TO^p = In^^-n^sl, where = dwp'^.iu) is 
the on-site particle density on each sublattice. The pos- 
sible phases of the Anderson-Hubbard model can then be 
classified as follows: The systems is a 

i) paramagnetic metal if iVf ""(O) ^ and mf^^ = 0, 

ii) AF metal if 7Vf°'"(0) ^ and mfp™ ^ 0, 



in) AF insulator if iVf °™(0) = and 



'AF 



but 



TVf °™(w) ^ for some w ^ 0, and 

iv) paramagnetic Anderson-Mott insulator if Nf''°™{uj) = 
for all uj. 

The ground state phase diagram of the Anderson- 
Hubbard model ll]) obtained by this classification is 
shown in Fig. [T] Depending on whether the interaction 
U is weak or strong the response of the system to disor- 
der is found to be very different. In particular, at strong 
interactions, U/W > 1, there exist only two phases, an 
AF insulating phase at weak disorder, A/W < 2.5, and 
a paramagnetic Anderson-Mott insulator at strong dis- 
order, A/W > 2.5. The transition between these two 
phases is continuous. Namely, the local DOS and the 
staggered magnetization both decrease gradually as the 
disorder A increases and vanish at their mutual boundary 
(lower panel of Fig. [2]). By contrast, the phase diagram 
for weak interactions, U /W < 1, has a much richer struc- 
ture (Fig. [1]). In particular, for weak disorder a param- 
agnetic metallic phase is stable. It is separated from the 
AF insulating phase at large [/ by a narrow region of AF 



FIG. 2: Upper panel: Staggered magnetization m|°p™ as 
a function of interaction U. Lower panel: niAp, oi — 
geom/arith, as a function of disorder A. Inset: Polarization 
Paf a function of disorder. Dashed lines present results 
obtained by arithmetic averaging. 

metallic phase. 

To better understand the nature of the AF phases in 
the phase diagram we take a look at the staggered magne- 
tization TO^p . The dependence of m^p"^ on U is shown in 
the upper panel of Fig. [2] for several values of the disorder 
A. In contrast to the non-disordered case a finite inter- 
action strength U > C/c(A) is needed to stabilize the AF 
long-range order when disorder is present. The staggered 
magnetization saturates at large U for both averages; 
the maximal values depend on the disorder strength. In 
the lower panel of Fig. [2] the dependence of m^p on the 
disorder A is shown for different interactions U. Only 
for small U do the two averages yield approximately the 
same results. Another useful quantity is the polariza- 
tion P^p = mlp/J«, where = E.. P"A^)d^/^ 
is the total spectral weight of p", (w). It allows one to 
investigate the contribution of the point-like spectrum 
of the Anderson localized states to the magnetization. 
This provides important information about the spectrum 
since with increasing disorder more and more one-particle 
states of the many-body system are transferred from the 
continuous to the point-like spectrum. For weak inter- 
actions {U = 0.5) the decrease of the polarization with 
increasing disorder A obtained with geometric or arith- 
metic averaging is the same (see inset in Fig. [2|) . Since 
arithmetic averaging does not treat states from the point- 
like spectrum correctly, the decrease of m^p (which is 
also the same for the two averages, see lower panel of 
Fig. m must be attributed to disorder effects involving 
only the continuous spectrum. At larger U the polariza- 
tion is constant up to the transition from the AF insu- 
lator to the paramagnetic Anderson-Mott insulator. In 
the latter phase the polarization is undefined, because 
the continuous spectrum does not contribute to /^™™. 

The AF metallic phase is long-range ordered, but there 
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FIG. 3: Typical local DOS as a function of disorder A for 
interaction [7 = 0.5 (upper panel) and IJ = 1.0 (lower panel). 
Solid and dashed lines represent opposite spin directions. 



is no gap since the disorder leads to a redistribution of 
spectral weight. In Fig.[3]the local DOS in the vicinity of 
the transitions between the paramagnetic metal, the AF 
insulator and the AF metal at [/ = 0.5 (upper panel), and 
the transitions between the AF insulator, the AF metal 
and back into the AF insulator at ?7 = 1.0 (lower panel) 
are shown. The paramagnetic metal, where p'^^iio) = 
p'^g{iu), is seen to be stable only for weak interactions. 

In the absence of disorder the AF insulating phase has 
a small ("Slater") gap at U/W < 1 and a large ("Heisen- 
berg") gap at U/W > 1. These limits can be described 
by perturbation expansions in U and 1/U around the 
symmetry broken state of the Hubbard and the corre- 
sponding Heisenberg model, respectively. In agreement 
with earlier studies [40| our results for ttiaf (upper panel 
of Fig. [2) show that there is no sharp transition between 
these limits, even when disorder is present. This may be 
attributed to the fact that both limits are described by 
the same order parameter. However, the phase diagram 
(Fig. [T]) shows that the two limits can be distinguished 
by their overall response to disorder. Namely, the reen- 
trance of the AF metallic phase at A/W > 1 occurs only 
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within the Slater AF insulating phase. 

The magnetic structure of the Anderson-Mott insula- 
tor cannot be determined by the method used here since 
it describes only the continuous part of the spectra and 
not the point spectrum. However, only the paramagnetic 
solution should be expected to be stable because the ki- 
netic exchange interaction responsible for the formation 
of the AF metal is suppressed by the disorder. This does 
not exclude the possibility of Griffiths phase-like AF do- 
mains 
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It is interesting to note that even the DMFT with an 
arithmetic average finds a disordered AF metal [32. 34 1. 
However, the arithmetically averaged local DOS incor- 
rectly predicts both the paramagnetic metal and the AF 
metal to remain stable for arbitrarily strong disorder. 
Only a computational method which is sensitive to An- 
derson localization, such as the DMFT with geometri- 
cally averaged local DOS employed here, is able to de- 
tect the suppression of the metallic phase for A/W > 1.5 
and the appearance of the paramagnetic Andcrson-Mott 
insulator at large disorder A already on the one-particle 
level. 

In conclusion, we computed the ground state phase 
diagram of the Anderson-Hubbard model at half filling 
within a non-perturbative approach which can treat in- 
teractions and disorder of arbitrary strength and is sen- 
sitive to Anderson localization on the onc-particlc level. 
For low disorder and weak interactions paramagnetic and 
antiferromagnetic metallic phases become stable, with a 
reentrant behavior of the latter phase. Slater and Hciscn- 
berg antiferromagnets can be distinguished by their very 
different response with respect to disorder. Experiments 
with cold fermionic atoms loaded into optical lattices will 
be able to test these predictions and check the accuracy 
of the theoretical approach employed here. 
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